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Метод динамического туннелирования 
изображений пятен лазерных пучков 

для улучшения качества прогнозирования 
их характеристик 


В статье рассматривается проблема обработки и прогнозирования динамических изображений в реальном 
масштабе времени. Предлагаемый подход формирования эталонного туннеля предполагает метод обработки 
информации, при реализации которого вычисляются центры пятен лазерных пучков с повышенной 
точностью, а также изображения трассы классифицируются с целью улучшения результатов дальнейшего 
прогнозирования их поведения. 


Введение 


На данный момент существующие магистральные системы передачи испытывают 
необходимость в существенном расширении канала. До сих пор данная потребность 
решалась за счет применения оптоволоконного кабеля, однако применение его на даль- 
них дистанциях сопряжено с многими технологическими сложностями, которые в той 
или иной мере приводят к увеличению стоимости строительства и эксплуатации подоб- 
ной системы. Одной из самых перспективных технологий как на больших, так и на 
малых дистанциях являются атмосферные системы передачи, позволяющие передавать 
значительные объемы информации посредством лазерного луча. Особенно перспектив- 
ными в данном направлении выглядят системы, позволяющие передавать информа- 
цию со спутника на Землю [1]. 

Для обеспечения спутниковой связи с помощью лазера необходимо постоянное 
позиционирование принимающего и передающего устройств с высокой точностью. При 
этом величина ошибки позиционирования не должна превышать тысячной доли градуса, 
искажения приёма будут настолько высоки, что связь не состоится. 

Поскольку под влиянием движения спутника, а также искажения, вызванного тур- 
булентностью воздушных масс, геометрические параметры и положение пятна лазерного 
пучка постоянно изменяются, необходимо постоянно в режиме реального времени про- 
гнозировать его положение. 

Цель данной работы — решить задачу эффективного прогнозирования определения 
центра пятенного объекта с помощью максимального использования его информацион- 
ных признаков, а также классифицировать кадры последовательности изображений пятен 
лазерных пучков (ПЛП), с целью фильтрации лазерной трассы от сильно искаженных дей- 
ствием помех изображений, тем самым сформировать туннель эталонных изображений. 

В данной работе предлагается определять центры ПЛП с помощью контурной 
ленты, образующейся вследствие использования градиентных масок или других методов 
определения границ на нечетких изображениях. После чего необходимо проводить 
формирование эталонной трассы с помощью классификации, основанной на формирова- 
нии динамических границ последовательностей изображений. 
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Контурная лента — область, образующаяся из элементов изображения, которая 
позволяет отделить объект от фона. Большинство методов выделения контурных призна- 
ков основывается на том, что разным объектам изображения отвечают области с относи- 
тельно одинаковыми значениями яркости. На границах же яркость заметно меняется. 

В качестве примера в работе используется алгоритм фильтрации по алгоритму 
Робертса. После наложения градиентной маски на изображение определяются внешние 
и внутренние границы объекта. Внутренняя граница задаёт базовый центр тяжести 
объекта, в то время как область между внешней и внутренней границей используется 
для его коррекции. 

Для нахождения центра тяжести контура необходимо найти площадь и статические 
моменты инерции контура, которые рассчитываются по формулам, изложенным в ра- 
боте [2]. Следующим шагом является вычисление площади контура. Поскольку контур 
можно рассматривать как многоугольник, то площадь контура, заданная своими верши- 
нами в порядке обхода, определяется в виде замкнутой кривой без самопересечений. 

Коррекция центра объекта происходит следующим образом. 

Используются только точки, лежащие в области между внешней и внутренней 
границей. Сначала находим приведенное максимальное значение яркости в контурной 
области. Каждой точке в данной области х(1‚)), где 1, / — соответственно абсцисса и ор- 


дината точки, назначается своя пара весовых коэффициентов влияния, К; и К’,. Далее 


определяются параметры коррекции и находятся координаты скорректированного 
центра, что также описано в работе [2]. 


Формирование границ туннеля 


Изменения положения пятна в трассе можно описать случайным процессом. 
Предположим, что определенная последовательность кадров трассы содержит лишь уни- 
модальные значения координат центров пятен лазерных пучков [3]. Таким образом, 
можно утверждать, что появление определенного значения координат ПЛИ можно опи- 
сать нормальным законом распределения с математическими ожиданиями 11, (1), м, (Г) и 


среднеквадратическими отклонениями (СКО) о’ (1), 5. (Е). Тогда данные значения мо- 


гут быть найдены по соответствующим формулам: 


Е (1) 


Ур, Хон, 
о,(®= пы = ы . (2) 


В реальных условиях центр ЛИ в смежных кадрах за счет воздействия окружаю- 
щей среды или других факторов может быть несколько смещенным, поэтому в дан- 
ном случае одного распределения становится недостаточно и его параметры нужно 
корректировать. 

При формировании и обновлении координат границ туннеля [4], [5] удобно исполь- 
зовать смоделированные одним нормальным распределением средние значения ко- 
ординат центров лазерных пучков. Тем не менее, если скорость изменения координат 
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центра пятна не нулевая и туннель имеет узкие границы, существует вероятность по- 
падания координат центров, важных для прогнозирования пятен, во внетуннельную 
зону, что приводит к увеличению погрешности прогнозирования. Чтобы этого избежать, 
необходимо корректировать границы туннеля в соответствии с изменениями координат 
центров лазерных пучков, при этом, исключая кадры с существенно искаженными цен- 
трами, не принадлежащими к туннельной зоне. Предлагается выделять смещения центров 
ПИЛИ с помощью туннеля и корректировать туннель в случае плавных смещений. При 
формировании начальных границ туннеля целесообразно проверить следующие условия: 


о (1-Р)> о. (3) 


В случае выполнения условия (3) необходимо дополнительно проверить условие: 


ах 

—|> О, 

Ру (4) 
где Р е [0,95-—0,99] — доверительная вероятность, С -— пороговое значение 


смещения. 

Невыполнение условия (3) свидетельствует о том, что данная координата отвечает 
модели с накопительными параметрами, а выполнение условия (4) указывает на наличие 
существенного смещения центра. При формировании границ туннеля в случае, если 
условие (4) выполняется, предлагается учитывать параметры данной точки с наимень- 
шим весовым коэффициентом. 

Достаточную информацию для этого обеспечивает детектор выхода за границы 
туннеля, следящий за попаданием координат центра лазерного пучка в туннель и филь- 
трующий одиночные помехи. При наличии такой информации обновление границ тун- 
неля происходит согласно формулам: 


Хх, = (1-А)х М 
У = (-Юу, +, 


е К (х, в [0х -м/2 [ха +%/2; +5] (у, ®[-55; у -и/2 [у +%/2;+0]), 
_ Кох, в [ха и/ 2; +%/2] бу, в [у -и/2; у. +/2], 


ОА АТ, (5) 
где х, — координата средней линии туннеля в 1-й момент времени, /1,, — матема- 


тическое ожидание х, ‚ и’ — ширина туннеля, Ах, К› — параметры, определяющие ско- 
рость обновления границ туннеля. Поскольку невозможно создать идеальный детектор 
границ туннеля, и с учетом возможности резкого изменения положения центра ЛП, 
не рекомендуется А; присваивать нулевое значение, т.к. при этом кадры с внетун- 
нельными центрами никогда не будут учтены. Это может привести к полной потере 
входной информации при импульсном изменении центра положения пучка или к необ- 
ходимости ручной коррекции положения средней линии туннеля. При больших значе- 
ниях А2 влияние текущего значения положения центра лазерного пятна будет настолько 
весомым, что пятна с координатами, близкими к значениям границ туннеля, и соответ- 
ственно, в большинстве случаев с большим значением шумов, будут также существенно 
влиять на среднюю линию туннеля. Рекомендованные экспериментально подобран- 
ные значения: А/= 0,09; А2= 0,35. 
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Формирование границ тунеля 


Рисунок 1 — Формирование границ туннеля. А — изменение координаты центра ЛИ 
в трассе; В — динамическая граница туннеля; С — статическая граница туннеля 


На рис. | показано динамическое изменение границ туннеля в соответствии с из- 
менением положения координат центра лазерного пучка во времени. Из рисунка 
видно, что количество точек, попадающих в туннельную зону, значительно больше 
в случае формирования динамических границ, чем статических. Данная методика дает 
возможность классифицировать изображения трассы на группы изображений, коорди- 
наты центров которых вошли в туннельную или во внетуннельную зоны, принимая 
во внимание информативные признаки обеих групп. Таким образом, при прогнози- 
ровании следующего кадра учитываются не только «эталонные» изображения и со- 
ответствующие им математические ожидания, а и возможность наложения шума на 
изображение. 


Результаты экспериментальных исследований 


Обработав значительное количество изображений (15 лазерных трасс по 1000 изо- 
бражений в каждой), можно отметить, что данный метод улучшает результаты алгорит- 
мов прогнозирования поведения изображений на 15-— 20 %, за счет использования 
дополнительно информации яркостей точек, находящихся на границе между фоном 
и объектом, и устранения сильно искаженных действием помех для кадров изображений. 

Результаты обработки трассы № 10 с разными значениями ширины туннеля 
приведены в табл. 1. 

В табл. | введены обозначения: \\ — ширина туннеля; Н„-— количество попа- 
даний в туннель (координаты центров объекта получены разработанным алгоритмом 
определения центров с повышенной точностью); Н.— количество попаданий в туннель 
(координаты центров объекта получены классическим методом); Р‚, Ре - процент 
попаданий в туннель; К - разность процентного попадания в туннель разработанного 
и классического методов. 

Существует большое количество способов прогнозирования [6] величины сле- 
дующего шага некоторого временного ряда, отличающихся по количеству необходимых 
вычислений, погрешностью прогнозирования и количеством шагов, на которые может 
прогнозировать данный алгоритм. Однако в случае работы в реальном времени на пер- 
вое место выходят скорость обработки и аппаратные затраты на выполнение алгоритма. 
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Таблица 1 — Зависимость количества попаданий в туннельную зону 
Разработанный и классический методы 


4К 


\ Н, Р.(%) | аниам |н. Р.(%) | ануа\ | в(%) 
0,2 491 24,7 219 472 23.8 228 1,0 
0,3 750 т 259 702 35,3 230 2,4 
0,4 992 49,9 242 922 46,4 220 3,5 
0,5 1202 60,5 210 1124 56,6 202 3,9 
0,6 1410 71,0 208 1304 65,6 180 5,3 
0,7 1553 78,2 143 1490 75,0 186 3,2 
0,8 1694 85,3 141 1650 83,0 160 5) 
0,9 1788 90,0 94 1750 88,1 100 1,9 

1 1860 93,6 72 1809 91,0 59 2,6 
т 1909 96,1 49 1868 94.0 59 2.1 
1,2 1924 96,8 15 1906 95,9 38 0,9 
ПЕ 1944 97,8 20 1925 96,9 19 1,0 
1,4 1969 99,1 25 1941 97,7 16 1,4 
1,5 1980 99,6 п 1952 98,2 п 1,4 
1,6 1983 99,8 3 1964 98,8 12 1,0 
1,7 1984 99,8 1 1975 99.4 п 0,5 
1,8 1987 100,0 3 1983 99,8 8 0,2 
1,9 : з к 1984 99.8 1 0,2 

2 ь ь : 1984 99,8 0 0,2 
р я : Е 1985 99,9 1 0,1 
2,2 - - и 1985 99,9 0 0,1 
2,3 я Е : 1986 99,9 1 0,1 
2,4 я - Е 1987 | 100,0 1 0 


В табл. 2 представлены данные о средней ошибке прогнозирования некоторых 
методов прогнозирования для известного метода предобработки и метода ПКЛ без 
туннелирования (бхср — средняя относительная погрешность метода, рассчитывалась 


между действительным и прогнозируемым значением). 


Таблица 2 — Средняя ошибка прогнозирования значения координаты Х изображения 
центра пятна лазерного пучка без туннелирования 


Название метода 


Относительная погрешность 
при прогнозировании на 1 шаг, бхср,% 


Известный метод 


Метод предобработки 


предобработки «контурной ленты» 

Числовые модели[7] 
1)Экспоненциальное сглаживание, 

) ц 1,05 0,93% 
аддитивная модель 
2) Экспоненциальное 

) ы 1,05 0,93 
сглаживание, модель Винтера 
3) Авторегресионная модель [8] 2,47 1,87 


Нейронные сети [9] 


структуры МР5-5-4 


4) Радиально базисная функция 0,78 0,67 
РБФ 55[10] 

5) Линейная сеть 55 [11] 0,76 0,62 
6) Многослойный персептрон 0,76 0,66 
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Данные алгоритмы были промоделированы в программе ЗТАТИЗТТСА 6.0 [12]. 
Как видно из табл. 2 для одних и тех же методов, метод «контурной ленты» пока- 
зывает результаты на 15 -— 20 % лучше, чем известный метод. Данные ошибки прогно- 
зирования при использовании алгоритма туннелирования на примере координаты Х 
представлены в табл. 3. 


Таблица 3 — Средняя ошибка прогнозирования значения координаты Х центра 
пятна лазерного пучка с туннелированием 


Относительная погрешность при прогнозировании 
на 1 шаг, бур, 
Название метода ы 
Известный метод Метод предобработки 
предобработки «контурной ленты» 
Числовые модели 
1)Экспоненциальное 
сглаживание, аддитивная модель 0,59 0,87 
2) Экспоненциальное 
сглаживание, модель Винтера 1,00 0,87 
3) Авторегресионная модель 2,05 1,46 
Нейронные сети 
4) Радиально базисная функция 
РБФ $5 0,75 0,64 
5) Линейная сеть $5 0.74 0.62 
6) Многослойный персептрон 0,74 0,63 
структуры МР5-5-4 


Из табл. 3 видно, что метод «контурной ленты» является лучшим с точки зрения 
предварительной обработки для прогнозирования центра лазерного пучка. Также можно 
отметить, что погрешность дхсь меньше при использовании алгоритма туннелирования и 
нейронных сетей, поэтому они являются более предпочтительными для применения по 
сравнению с численными методами. 


Выводы 


К основным преимуществам данного метода по сравнению с классическими ме- 
тодами [13], [14] определения центра фигуры [15] можно отнести следующие: 

1. Точность определения центра увеличена за счет учета особенностей кон- 
турной линии. 

2. Алгоритм реагирует на наименьшие изменения границы объекта динамичес- 
кого изображения и соответствующим образом влияет на определение его центра. 

3. Формирование границ туннеля не требует значительных затрат времени. 

4. Диапазон допустимых входных значений ограничен лишь возможностями вы- 
числительной системы (алгоритм корректно работает с резкими изменениями поло- 
жения центра ЛП). 

5. Эксперименты показывают, что предложенный метод предварительной обра- 
ботки изображений на основе «контурной ленты» является более предпочтительным по 
сравнению с применением известных методов предварительной обработки изобра- 
жений для прогнозирования положения лазерного пучка. 
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